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in ' Abstract. We introduce a new fractional oscillator process which can be 

t*H^ obtained as solution of a stochastic differential equation with two fractional 

■^r 1 orders. Basic properties such as fractal dimension and short range dependence 

' of the process are studied by considering the asymptotic properties of its co- 

\I ' variance function. The fluctuation-dissipation relation of the process is inves- 

vH ^ tigated. The fractional oscillator process can be regarded as one-dimensional 

fractional Euclidean Klein-Gordon field, which can be obtained by applying 
the Parisi-Wu stochastic quantization method to a nonlocal Euclidean action. 
The Casimir energy associated with the fractional field at positive temperature 
is calculated by using the zeta function regularization technique. 



1. Introduction 



Work on extending stochastic processes characterized by a single index to corre- 
sponding processes with two indices or processes with variable index has attracted 
considerable interest lately. For example, fractional Brownian motion parametrized 
^ ' by a constant Hurst index H has recently been generalized to the bifractional Brow- 

^^D , nian motion [HE] characterized by two indices, with fractional Brownian motion as 

^^ ' a special case. A process with two indices provides a more flexible model as far as 

applications are concerned. Fractional Brownian motion has also been extended to 
multifractional Brownian motion with a variable index H{t) [2111]. Multifractional 
Brownian motion allows one to model phenomena with variable memory and fractal 
dimension that changes with time or position. Another example is the process of 
f^ I Cauchy class [51, which has been extended to the process with two- indices known 

as generalized Cauchy process [B] [7] . In contrast to many other single-index pro- 
cesses (for example, fractional Brownian motion), generalized Cauchy process has 
the advantage that its two indices separately characterize the fractal dimension or 
self-similar property and the long range dependence. The aim of this paper is to 
introduce a new type of fractional Gaussian process indexed by two parameters so 
that its short range dependence property and fractal dimension can be separately 
characterized. Such a process can be regarded as a generalization of fractional 
oscillator process or fractional Ornstein-Uhlenbeck process with single index [S] [S] . 
This paper is organized as follows. Section 2 introduces the generalized fractional 
oscillator process as solution of a stochastic differential equation with two fractional 
orders. Despite the covariance of this process does not have a closed analytic form, 
its basic properties can be studied by considering the asymptotic properties of 
its covariance. In particular, the fractal dimension and short range dependence are 
studied. The fluctuation-dissipation relation are discussed in section 4. In section 5, 
we consider the fractional oscillator process as one-dimensional Euclidean fractional 
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scalar field. Stochastic quantization of the field at zero and finite temperature 
is carried out. We proceed to calculate the Casimir energy associated with the 
fractional field at finite temperature by employing the zeta function regularization 
technique. In the last section, we briefly discuss possible generalizations of the 
results obtained. 

2. Fractional Oscillator Process with two indices 

In this section we define the fractional oscillator process with two indices. Re- 
call that an ordinary oscillator process X{t) can be obtained as a solution to the 
Langevin equation 

(2.1) iDt + X)X{t)^7^it), 

where 7y(i) = dB{t)/dt is the standard white noise with 

(r;(i))=0, {Ti{t)r^{t'))^6{t-t'), 

B{t) is the standard Brownian motion, and A is a positive constant. Using Fourier 
transform, the solution of (|2.ip can be written as 

X{t) = -= / - — ^dw. 



where ??(w) is the Fourier transform of the standard white noise ri{t): 

flic.) ^ ^= / e-^'-dB{t). 

In the literature, X(t) is known as the oscillator process or the Ornstein-Uhlenbeck 
process. It is a centered stationary Gaussian process with covariance function 

(2.2) (x(s)X(s + i)) = — y -^-^dc^- 



2A 

One can also regard the oscillator process X(t) as one-dimensional Euclidean scalar 
Klein-Gordon field with mass A, and propagator given by the spectral density 

1 1 



5(6.) 



27r w2 + A2 ■ 

Since fractal dynamics [10[ lll[ 112] have increasingly played an important role in 
various transport phenomena in complex media, it would be interesting to inves- 
tigate various possible generalizations of X(t) to its fractional counterpart. One 
direct way is to replace the differential operator Dt in (12. 1|) by the fractional dif- 
ferential operator a^t to obtain the following type I fractional Langevin equation 

[Tana [151 [16]: 

(,Df + A)X„a(t) = r;(i), 
where the fractional derivative a^t is defined as [mfTSlfTO] : 

( aD-^!) (t) = —^ (4-) I 7 ^1"^ ^, du, n-l<a<n. 

When a — 0, qD" is known as the Riemann-Liouville fractional derivative; and 
for a — — cx), -oo-Df is called the Weyl fractional derivative. Another possible 



FRACTIONAL OSCILLATOR PROCESS WITH TWO INDICES 3 

generalization is to fractionalize the operator {aDt + A) to obtain the foUowing type 
II fractional Langevin equation [H [9] : 

UA + A)^Xi,^(t)=r7(i), 7>0. 

Recently we have combined the investigation on these two type of processes and 
study the more general case pUj : 

(2.3) {aD^ + XyX^^^{t) = 7j{t). 

If Weyl fractional derivative is used in (|2.3p . then for aj > 1/2, Xa^^{t) turns out 
to be a centered stationary Gaussian process with a representation 

1 f e'"*77(^) ^ 



and covariance function 

(2.4) 

(X„,,(.)X„,,(. + t))^-y ,(^^)„^„..rf- = ^y (M- + 2AMcos^ + A^)^'"- 

The properties of the process Xq ,y(t) have been studied in [20]. However, the above 
generalization contains an unsatisfactory aspect, namely its spectral density 

(2.5) S(uj) = —- ^ -^ 

^ ' ^ ' 2^(|cj|2" + 2A|w|cos2^ + A2)^ 

has a complicated form. When a — 1, the spectral density simplifies to 

Siu;) ^ ^ 



27r(|a;|2 + A2)T- 
This signifies that Xi^-y(t) has another representation 



^/2^i (|c.|2 + A2)^ 
and therefore is a solution of the following fractional stochastic differential equation 

(-A + A2)*Xi,^(t)=ry(i), 

where A = d?/dt^ is the one-dimensional Laplacian operator. Inspired by this, 
we define a new type of stochastic process Xa^j{t) indexed by two parameters a, 7 
with a E (0, 1], 7 > 0, as a solution to the following fractional stochastic differential 
equation: 

[(-A)" + A2]*X„,^(t)=r7(t), 
or equivalently, the equation 

(2.6) (D^ + A^)*X„,^(t)=rKi), 

where Dj " is the Riesz derivative defined by [ITl [21] 

D^/:=(-A)"/:=F-i(|c.|2"/H), 
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C(=0.4. p=1 .5, ?.=2 



;=0.6, P=1.5.X=2 





Figure 1: Some simulations of Xa.-iit) for different values of a and 7. 



with / := F{f) the Fourier transform of /. It follows easily that when 07 > 1/2, 
the solution to (|2.6p is given by 



(2.7) 



Xa^^yi) 



1 



2^1 (\u^\2^ + x^y 



-doj. 



We also call Xa^'^{t) a fractional oscillator process. Its covariance function Ca,'i{t) — 

{Xa^^{s+t)Xa,'i{s)) is 

(2.8) c^,^{t) ^^ r 



jdoj. 



(|w|2"+A2)^ ' 

and the spectral density has the required simpler form 

1 1 



(2.9) 



Sa,jiuj) 



2tt {\uj\^" + x^y 

compared to (|2.5p . Some simulations of the process Xa,'y{t) are given in Figure 1. 

Here we would also like to remark that when 127 < 1/2, Xa^^yit) only exists as a 
generalized stochastic process, in the sense of generalized functions. Namely, for a 
Schwarz class test function f{t), 



{X^,^{t)J{t)) 



(k 



2a 



2A2 



A2) 



— du), 



FRACTIONAL OSCILLATOR PROCESS WITH TWO INDICES 5 

where /{lo) is the Fourier transform of /. 

The two fractional generahzations Xa^j{t) and Xa.-yit) of ordinary osciUator pro- 
cess do not give the same solutions since they give rise to different spectral densities. 
However, if we want to regard fractional oscillator process as one-dimensional Eu- 
clidean Klein-Gordon field, (|2.6p is one-dimensional Klein-Gordon equation with 
two fractional orders, and the covariance function (|2.8p becomes the propagator of 
the corresponding fractional Klein-Gordon field (f)a,^{t). We shall show in section 
5 that (j)a,-y{t) can be obtained by applying Parisi-Wu stochastic quantization [22] 
involving a nonlocal Euclidean action. 



3. Asymptotic behaviors of the covariance function C'a^-yit) and 

SAMPLE PATH PROPERTIES OF Xa,^{t) 

3.1. Asymptotic behaviors of the covariance function Ca^jit). When a — 1, 
7 > 0, the covariance function Ci^^(i) (|2.8p has the following closed form [8]: 

91/2-7 /|i|\'''"l/2 

(3.1) Ci.^(i) = -^=^ (^Uj X^_i/2(A|t|), 

where K^{z) is the modified Bessel function of second kind or the MacDonald func- 
tion. However, the covariance function Ca^jlt) in general does not exist in closed 
analytic form. In fact, since the spectral density (|2.9p has the same functional 
form as the characteristic function of the generalized Linnik distribution ^3" , the 
covariance function Ca^-yit) (|2.8p has the same functional form as the Linnik prob- 
ability density function, whose analytic properties have been studied in [24l [23] . In 
particular, one can obtain the following integral representation for the covariance 
function Ca^^(t): 

(3.2) C^^^{t) = -Im 



'^ Jo (A2 -he-""u2")''' 

It turns out that the analytic properties of Linnik probability density function 
depend on the arithmetic nature of the parameters a and 7; and the conditions im- 
posed on a and 7 are rather complicated and are not of practical interest. Therefore, 
we shall use different methods to study the asymptotic behaviors of the covariance 
function Ca.j{t) that are more suited for various applications. 

For the properties of XQ,_^(t) that we are interested such as its fractal dimension, 
long or short range dependence, it suffices for us to know the leading behavior of 
the variance of the associated increment process 

(3.3) a^t) = ([Xa,^(s + t) - A„,^(s)]') = 2Ca,^(0) - 2C„,^(i) 

for t ^ and the leading behavior of Ca,-y{t) for i — > 00. 

First, we examine the behavior of (T^(i) when t — > 0. We impose the restriction 
aj > 1/2 so that Xa,j{t) has finite variance and p.3p is well-defined. Under this 
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restriction, the variance Ca,-y{0) is given by (#3.251, no. 11, Ref. [25] ) 

2T:a r(7) 

2ar(7)r(l-7+^)sin(.(7-5k))' 

where the identity —zT{z)T{~z) — 7r/sin(7rz) has been used. Our result is in 
agreement with that of Erdogan and Ostrovskii [53] . To obtain the leading behavior 
of 

(3.5) a'{t) = 2C„,,(0) - 2C.,,(i) = |^°" l^S^^du;, 

we consider the cases 1/2 < a'f < 3/2, aj > 3/2 and aj = 3/2 separately. 
Case I. When 1/2 < aj < 3/2, ((S31) is equal to 

(3.6) a^{t) — / ,,, \ ( d^ 
4|^|2a7-i /-oo sin2(a;/2) , 



TT 7o (C^2a^.^2|i|2a)7 



^'*''"'' ' '"a;-2a7sin2(|)dc. + o(|i|2"^-i) 



TT 







'' +o(|t|2"T-i) as i 



cos(7ra7)r(2a7) 

We have used #3.823 of '25] in the last step. Eq. ([3?6)) shows that when 1/2 < 
aj < 3/2 and i — > 0, the leading term of a'^{t) is of order |i|2"T~i. If we replace 
a J by H + 1/2, then (j3.6p becomes 

U|2H 

sm(7r//)l (2ii + 1) 

which shows that the short time asymptotic behavior of (T'^{t) is characterized by 

the index H = 07 — 1/2. 

Case 11. For aj > 3/2, using 1 - cos(w|i|) = ^cj^jip + 0{\t\*) as t ^ 0, we have 

This shows that after crossing the point a'j ~ 3/2, the leading behavior of o'^(t) is 
of order jtp, which does not depend on the parameters a and 7. 
Case III. The borderline case aj = 3/2 is more complicated. First, we find as in 
Case I that 
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The integral 

3.8) / ^^ \' ' , dw 

does not converge when t ~* because of the singularity at the origin of the 
integrand iu~^°'^ sin^ {^/'^)- Since sin 2; r-^ z when 2; — > 0, we write p.Sp as a sum of 
two terms Ai{t) and A2(t), where 



•^■'"^iV 



duj. 



+ A2|f|2")7 

reflects the divergence of p.Sp when i — > 0; and 

/•I sin2(co/2) - (co/2)2 , /-"^ sin2(w/2) 

^' Jo (tj2"+A2|i|2")7 J^ (^2a_,_^2|i|2a)7 

carry the finite part. By making a change of variable v = uP''^ or equivalently 
ui — w'*'/'^, we find that 

1 f^ v'^-^dv 



^'W^s^i 



8q7o (W + A2|i|2")7- 

To find the asymptotic behavior of Ai{t) as i ^ 0, we split it again as the sum of 
Az{t) and A4^{t), where 

1 /-^ dv 111 

8^ io (w + A2|i|2") ^ 4 °^ ^ ~ 4^ " 
give the divergence part; and A4^{t) :— Ai{t) — A^{t) gives a finite limit: 
(3.10) 



(3.9) A3(i) = — / , , ^,,,, . ^ - log - - — log A + 0(1) 



,,, 1 fl w'^-l-(l' + A2k|2")7-l 1 /-IT^^ „7-l_ (l+^,)7-l 

^' SaJo (« + A2|i|2")7 8a 70 (1 + ^^^ 

' ^ ' -d« + o(l) = — (^(l)-V(7)) + o(l). 



8a Jo (1 + '")'' 8a 

In the last equality, we have used #3.219 of [23], with ^p{z) = r'{z)/T{z) being the 
logarithmic derivative of the gamma function. The limit of A2{t) when t ^ is 
given by 

It can be computed explicitly by regularization method: 

(3.11) yl2(0) = lim^ | / 10^^+" sin^ (-) duj - - f cj-^+'duj 

,. (1 Tre r(l+e) 1 1 1 /,,,, 3 

Combining (|3.9p , (|3.10p and p. lip , we find that the leading behavior of a"^ (t) when 

aj — 3/2 is: 

(3.12) 

a\t) ^ ^ log 1 - ^ (i log A + ;! (V(7) - ^(1)) - ^(1) - I) + 0(1). 
n \t\ n [a 2a 2 J 
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Figure 2: These graphs show the small time behavior of cr^(t). Left: The func- 
tion — cos(7ra7)r(2Q7)(7^(f) is plotted as a function of t. The reference curve is y = 
|tp"^-\ Here 07 = 0.6. The graph shows that - cos(7ra7)r(2Q7)cr^(f) ~ \t\'^°'''-'^ 
when t ^ 0. Right; The function K^.ya^{t) is plotted as a function of t, where 



K^^^ = 27vaX^''' 



r(7) 



Ka,~,cy'^{t) ~ t^ when f ^ 0. 



The reference curve is y 



The graph shows that 
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Figure 3; The graph shows the small time behavior of a^it) when a'y — 3/2. The reference 
curve is t^ log(l/|f |). 



This shows that the leading behavior of (T^(i) at the borderline case 07 — 3/2 is of 
order \t\'^\og{l/\t\). 

The behavior of a'^ (t) at small t is illustrated graphically in Figure 2 and Figure 
3. We can further confirm our results by checking the small time behavior of (T^(t) 
for the case a — 1, with the covariance function Ci_^{i) given explicitly by (|3.ip . 
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By using #8.445, 8.446 and 8.485 of Ref. ^ one gets for v (^7L, 



IT Jy. (Z/2)^J-'- y. (Z/2)^J' + - 

2sin(^i.) ] f- j!r(j + l-v) f- j\T{j + l + v){' 



(3.13) KM = K^M - ^.-^ \ E MH^ - E 



and in the case v = ±77i, ra a nonnegative integer, 
(3.14) 

1 '^ (-l)J(TO-j - 1)! /Z\2i-™ 



^'>(-) =:j E 



2 ^ j\ 



+ (-ir+^E%^^^|iogf-^^(j + i)-iv'(j- + i+m)j, 

^j!(to + j)! I^ 2 2 2 J 

From ()3.13p and (I3.14p . one finds that the variance of Xi.^(t) is 



Ci,^(0) 



2A27-isin(7r(7-i)r(7)r(|-7)' 



and the leading behavior of (T^(i) as i ^ is given by: 

• If 1/2 < 7 < 3/2, 

^ ^ ^ ''^'^^22.-isin(.(7-i))r(7)r(7-i) 

cos(7r7)r(27) ' 

In the last step, we have used the identity T{2z) = {2"^'--'^ / ^/t:)T{z)T{z + (1/2)). 

• If 7 > 3/2, 

(^•^«) ^^W-4^^'*'^- 

. If 7 = 3/2, 

(3.17) a2(t)^^(^logl-logA + log2 + V(l) + i 

By putting a = 1 in the small time asymptotic formulas p.6p . (|3.7p and p.l2p 
of cr^(i), one obtains the formulas (|3.15p — (|3.17p . which confirm our results. 

Next, we study the asymptotic behavior of Ca^jit) for t -^ cx3. When < a < 1, 
we can make use of p.2p . Making a change of variable u i— > u/t and using the 
formula 

1 _^(_ipr(7 + j) , 



E 



(1 + z)7 ^^ j!r(7) 



z^ 
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Figure 4: This graph shows the large time behavior of Ca,-y{t). 



we can derive the following t — > oo asymptotic expression for Ca,-y{t) which is valid 
when a S (0, 1): 

(3.18) 



c„,T,(t) ^-imlr^ I 



e "du 



(A2 + e-*^"f^) 



■-Im It-' 



^u g tm:t±Ae-^^-j^^-2^-2, 



J=0 



j!r(7) 



f2aj 



7rr(7) 



J2 ^ r(7 + j)r(l + 2aj) sin {iraj) ^-(^"J+i). 



j=i 



The leading term is 

(3.19) Cc^At) - -A-2('^+i)r(l + 2a) sin (ttq) i-(2"+i). 



By letting A = 1, (|3.18p is in agreement with the analogous result given in Ref. 
[23] for Linnik distribution. Notice that when < a < 1, Ca.-f{t) is of polynomial 
decay with order ^-^a-i .^^j^gj-^ ^ ^ qq (■ggg pigurgg 4 and 5). 

The case a — 1 has to be considered separately. Using the relation (#8.451, no. 
6, ref. pF): 

4(^2 _ 1 



K,{z) 



1 



2z V 8z 

we find from the explicit formula (j3.ip for Ci^^(i) that 



(3.20) 



Ci,^(i) 



itri 



-A|t 



as Itl 



oo. 



One notice that at large time, Ci^^(i) decays exponentially, in contrast the large 
time behavior of Ca -yit), a G (0, 1) (|3.18p . which decays polynomially (see Figure 
5). 
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Figure 5: The large time behavior of C{t) when a — 0.6 and a — 1 respectively. Left: 



Here K — ^ ^ 

J.J.C1C .fVa,^ 7r(l + 2a)sin(ira) 

that K^_^Cc,-y{t) ~ f-2"-i 
the reference curve is y = e 



and the reference curve is y = t 



-(2q+1) 



The graph shows 



as t 



-^ oo. Right: Here N^{t) = {2\yT{'y)V-^ , A = 1 and 
The graph shows that iV.^(t)Ci,^(t) ~ e"^' for large t. 



From the above results, it appears that the small time asymptotic behavior of the 
covariance Ca^-yit) varies as |^|min{2a7-i,2}^ depending on both a and 7. However, 
if the index 7 is replaced by 7/a, then C„,^(i) - |i|nim{27-i,2} as i ^ 0. Thus, 
together with the large time asymptotic behavior Ca.-i{t) ~ |t|^^"^^ as t -^ 00, we 
have the result that the small and large time asymptotic behavior of the covariance 
of Xa^-y{t) are separately characterized by 7 and a. The physical implications of 
this result will be discussed in the subsequent sections. 



3.2. Locally Asymptotically Self-Similarity and Fractal Dimension. Recall 
that a stationary random process cannot be self-similar ^26] . It would be interesting 
to see whether Xa.-^it) satisfies a weaker self-similar property, namely self-similarity 
at very small time scales. First we introduce some definitions. A positive function 
/ is asymptotically homogeneous of order k at 00 if there exists a non-zero func- 
tion /°° such that, for almost every cj e K. and r > 0, /'^(w) — r^'^ J{ruo) has a 
limit f°^{u!) when r -^ cx). Clearly, f°°{uj) is homogeneous of order k, thus fixes 
the index k uniquely. One can easily verify that the spectral density Sa,j{uj) is 
asymptotically homogeneous of order 2a7 at 00, with S'^^(ijj) — cj^^"'^ / {2tt) . In 
addition, the spectral density satisfies the following property: there exist positive 
constants A, B e R such that Sa,^{Lu) < B\uj\~'^"'^, for almost all \lu\ > A. This is 
clearly true since A > implies 



Sa^-yiuj) 



1 



1 



-2q7 



27r(A2-|-w2a)7 



< 



2tt 



Using a result (Proposition 2 in Ref. [23), one concludes that if aj G (1/2,3/2), 
the fractional process Xa^^{t) is locally asymptotically self-similar (LASS) of order 
aj — 1/2, that is for w S M, 



lim 

6^0 + 



Xa,i{to + Su) — Xa,j{to) 



{^r,(")}, 
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Figure 6; The graphs shows that e "i'+2 [Xa.-f{et) — Xc,^(0)] approaches the fractional 



Brownian motion B _i{t) when e -^ 0. Here Zi{t) = e^ 
ei = 0.0001, £2 = 0.00005, J = 1, 2. 



-a7+o 



[Xa,-,(eit) - X^,^{0)] with 



with the convergence in the sense of distribution on the space of continuous paths on 
M. Note that this result is in agreement with (|3.6p which asserts that the covariance 
(J^it) ~ Cjtp"'''"^ as t -^ 0. The limit process or the tangent process X^ is self- 
similar of order a'f — \- It can be identified with fractional Brownian motion if 
a'y = H -\- 1/2, where H denotes the Hurst index of the fractional Brownian motion. 
Just like the case of ordinary oscillator process, which locally behaves like Brownian 
motion, likewise the fractional oscillator process Xa,j{t) has the same local behavior 
as fractional Brownian motion of order aj — 1/2. In fact, (|3.6p gives us 



lim 

r^0+ 



Xa,j{to + £u) — Xa^j{to) 



Xa,-,{to + ev) - Xa,y{to) 



= lim 

e— 



a'^ieu) + (T^(£u) — cr^(e(u — v)) 



2£2"7- 



(h 



|2a7— 1 



|2a7— 1 



|2q7~1'i 



2 cos(7ra7)r(2a7) 
which is the covariance function {Bh(u)Bh{v)) for the fractional Brownian motion 



Bniu) 



1 



- 1 



\H+h 



duj. 



if we identify H with a'y — 1/2 (see Figure 6). We also remark that when a/3 — 3/2, 
we find from (I3.12p that the fractional oscillator process Xa^-yit) fails to satisfy the 
LASS property. When wy exceeds 3/2, the process Xa^pit) becomes differentiable 
with variance 



lim 



^0,7(^0 + £) ^ ^0,7(^1 



n 2\ 



Q,7V''0j 



= lim 



.^(£)_A^-rr(^)r(7-|,) 



27ra 



r(7) 



which follows from p.7p . 

Another important concept in the study of the sample path properties of a 
stochastic process is the Holderian property. A function / : [a, 6] ^ E is Holderian 
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of order k G (0, 1] if 

l/W - /(s)| < K\t ~ sr, for all s, t e [a, b] 

for some constant K > 0. It is well-known that if Z{t) is a stationary process and 

cr2(t) = {[Z{t) - Z{0)f) satisfies 

a^t) <C|t|2-, 

then almost surely (a.s.) the sample path of Z{t) is Holderian of order k ~ e for 
all £ > [28l [29l |30] . Applying this concept to the fractional oscillator process 
^a.-yit), we find from p.6p . p.7p and (|3.12p that for any e > 0, the sample path of 
^a.-y{t) is Holderian of order 07 — 1/2 — e ii a-f < 3/2, and is Holderian of order 
l-£if a7>3/2. 

Now we want to consider the fractal dimension D of the graph for the fractional 
oscillator process Xa.-y{t). Since fractal dimension is a local concept, fractality is 
defined for infinitesimally small time scales. For a locally self-similar process, one 
may apply the following result to obtain its fractal dimension. A process which is 
LASS of order k > and its sample paths are a.s. k— e-Holderian for all e > 0, then 
the fractal dimension of its graph is a.s. equals to 2 — k [551[3T]. Applying this result 
to Xa^^{t) gives the fractal dimension D — ^ — aj for the graph of the fractional 
oscillator process when 1/2 < a'j < 3/2. For aj > 3/2, the fractal dimension 
of the graph of Xa,j{t) is equal to 1. In other words, the fractal dimension of the 
graph oi Xa^-yit) is max{l, | — 07}. Again, if we replace 7 by 7/a, then the fractal 
dimension becomes max {l, | — 7}, which depends solely on 7. 

3.3. Short Range Dependence Property. First we recall that the Ornstein- 
Uhlenbeck process or ordinary oscillator process (up to a multiplicative constant) 
is the only stationary Gaussian Markov process, thus rule out the possibility of the 
fractional oscillator processes Xa,'y{t) being Markovian. Now we want to find out 
the nature of memory possessed by Xa^-yit), whether it has long memory or long 
range-dependence (LRD), or short memory or short-range dependence (SRD). A 
stationary Gaussian process with covariance C{t) is said to be LRD if for some 
finite t, > 0, 

/•oo 

(3.21) / \C{t)\dT^ 00, 

otherwise it is SRD. From the results (|3.19p and (|3.20p . we find that the covari- 
ance of Xa,-y{t) behaves asymptotically as Ca,-f{t) ^ ^-(2"+!) if q: g (Q, 1) and as 
Ca,-y{t) ~ e~^^V~^ if a = 1, for t — > 00. This shows that the corresponding integral 
(|3.2ip is convergent and therefore the fractional oscillator process ^^.^(i) has SRD. 
It is interesting to note that when a € (0, 1) the asymptotic order of the covariance 
does not depend on the parameter 7. Therefore, the parameter a characterize the 
asymptotic order of the covariance Ca.-yit) as i — > 00. Combining with the remark 
given earlier, one notes that it is possible to separately characterize the fractal 
dimension and short range dependence oi Xa,-y{t) with two different indices. 

4. Fluctuation-Dissipation Relation 

One of the important theorems in statistical mechanics is the fluctuation-dissipation 
theorem |32] which relates the coefficient of the covariance of the external random 
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force in the Langevin equation with the frictional coefficient. It would be interest- 
ing to see whether there exists some kind of fluctuation-dissipation relation for the 
fractional process Xa^jit). For this purpose we re-express the fractional Langevin 
equation (|2.6p as 

(4.1) (D2" + A2")*X„,^(t) = r;(t), 

with the covariance of white noise ri{t) given by (77(^)77(5)) ~ 2BS{t — s), where 
S is a constant coefficient. If wc regard Xa^-y{t) as the velocity process, then the 
assumption of the thermalization of the fractional velocity process based on the 
fractional generalization of the classical principle of equipartition of energy [32] 
gives 

(4.2) ([Xo.,^{t)f) = ikTr\ 



where we have assumed the particle under consideration has unit mass; k is the 
Boltzmann constant and T denotes temperature. Under the condition 07 > 1/2, 
the process Xa^^{t) has finite variance given by (|3.4[) : 



(4.3) Ca^AO)— \./.l ^°^ A^-^"^. 



na r(7) 



From (14.21) and (14.31) one obtains 



(^■4) ^ = wXTTT^^^XY^'"^~'(^^)"^ = n(a,7)A2"'^-i(fcr)"^ 

^ \2a) ^ U 2a) 

with 

7rar(7) 



n(a,7) 



r(^)r(7-5k] 



can be regarded as the generalized fluctuation-dissipation relation for the 
fractional velocity process. When a = 1 and 7 = 1, (|4.4p reduces to the fluctuation- 
dissipation relation for the ordinary Ornstein-Uhlenbeck process Xi,i(i): 

B = XkT. 

We can show that for a short range dependent process such as Xa^-yit), the lead- 
ing term for the large time behavior of the variance of its mean-square displacement 
does not depend on the covariance of Xa^jit). For illustration, let us first consider 
the simple case with position process Ya.-y{t) linked to the velocity process by 

(4.5) X^,,it) = ^^^. 

If we assume that Ya^j{0) — 0, then the variance of the position process Ya^^{t) is 
given by 



(4.6) {[yan{i)f)^l f Cc,^i\si-S2\)ds2ds 
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By some calculus, we have 

pt ps-i pt pt 



(4.7) {[Ya,y{t)f\ ^ / Ca.-yisi- S2)ds2dsi+ j / Ca.^{s2- Si)dsids2 

^ ' Jo Jo Jo Jsi 

pi ps pt i*S2 

= / Ca.',{T)dTds + / / Ca,jis2 — Sl)dsidS2 

Jo Jo Jo Jo 



10 Jo Jo Jo 

= 2 / / Ca.,-y{T)dTds = 2 f {t - T)Ca.j{T)dT. 

Ja Jo Jo 

Since the integral L Ca^'y{T)dT is convergent, we find that in the long-time < 3> 1 
limit, 

r /"^^ 

t 



(4.8) 


([r„^^(i)]')~2 


/ C^rfir)dT 
Jo 


which 


is just ordinary diffusion with diffusion constant 




D= j C^,^iT)dT. 

Jo 


Since 






poo 


(4.9) 


2 / C„,^(T)dT = 2 


7rS'„,-,(0) = 2BX 



using (|4.4p . with {kT)'^^ replaced by kT in view of (|4.5p . one gets 

7rar(7) kT 

which reduces to the well-known Einstein relation 

for a = 7 = 1. This simplified example shows that the long time behavior of Ca.-f{T) 
does not show up in the leading term of the long time asymptotic expression of 
the variance ([yQ^^(i)] V Its effect only appears in the second leading term (see 

appendix). This is due to the fact that Xa,'y{t) is a short range process with its 
covariance C{t) ~ r"'', k > 1 for t — > oo. We remark that (|4.8p is consistent with 
the result obtained for the case with a > 0, 7 = 1 if the usual velocity-displacement 
relation (|4.5p is used [13]. 

Now we consider the fractional case with the velocity linked to the displacement 
by the following relation: 

(4.10) X^,^{t) = oD'tY^.^-.it), \<X<\. 

where ^Df denotes Riemann-Liouville fractional derivative of order x- If we further 

assume that oDf~''Ya ^U) = for j = 1 if x < 1 and j = 1, 2 if x > 1, then 

' t=o 
the position process is given by 



1 '■' 



(4.11) y„,^(t) = oi^x^^^it) ^Yr~) I (^ " u)^''X»M^)d 



lu. 
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One can show that (see appendix) 
(4.12) 



[r„^^(t)]') = 2i?A-2"7 



^Q|^^max{0,2x-2.2x-2a-l}i^g^ 



.(2x-i)r(x)2_ 

Note that the term in the bracket on the right hand side of (|4.12[) is just the variance 
of the Riemann-LiouviUe fractional Brownian motion indexed by x — 1/2 33J. By 
setting x = oil ^^^d substituting B from ()4.4I) . we get 

2\ 27rar(7) (fcr)"T f f°'^- 



(4.13) {[Yc.At)\ 



r(2k)r(7-5k) ^ 

^7V(a,7)('^^ 



(2^7 - l)r(a7)2 



a7 



{-) ' 



2a7-l 



where N{a,"f) is a constant term depending on a and 7, which can be normaUzed 
to become unity. Let P{x, t) denotes the probabihty distribution of the equihbrium 
displacement process (i.e. when t -^ 00), then the effective Fokker-Planck equation 
is 

lpi:,,t) = Dit)-^P{x,t), 
with the diffusion coefficient D{t) given by [13] 

i^(t)4|([r.„(t)f). 
For equilibrium state, \[Ya.-y{t)] ) is given by (I4.13p . so 



which can be regarded as the effective diffusion coefficient for fractional Brownian 
motion if we let aj — 1/2 = H, the Hurst index [34, 35, 36 . Our result differs from 
that of ref. [Ml EH [35] which has (fcT/A) instead of the {kT/XY"' given above. 
This is due to our use of A" in the Langevin equation (|4.ip and the fractional 
generalization of equipartition principle (j4.2p . 

From the discussion above, we see that the long-time dependence of the covari- 
ance of Xa^^(t) does not enter in the leading term of the variance of Ya^-y{t). It only 
appears as second leading term. Instead, the property of the leading term depends 
on the differential relationship between Xa^-^it) and Ya,-y{t). In the following, we 
return to the case where {r]{t)'q{t')) — 5{t — t'). 

5. Casimir Energy associated with Xa.'^[t) at Finite Temperature 

The ordinary oscillator process can be regarded as one-dimensional Euclidean 
scalar massive field as its spectral density (w^ -I- A^)""'^ is just the Euclidean prop- 
agator for a scalar field with mass A. By analogy, one can consider Xa^^{t) as a 
fractional Euclidean scalar massive field in one-dimension with propagator (|a;p" -|- 
A^)"'''. From this viewpoint, it will be interesting to find a convenient way to quan- 
tize the corresponding fractional quantum field (j)a^-y{t). This can be achieved by 
using the stochastic quantization of Parisi and Wu [52]. According to this quan- 
tization scheme, an additional auxiliary time t is introduced and the Euclidean 
quantum field <pa,-yit',T) is assumed to evolve in this auxiliary time according to 
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a stochastic differential equation of Langevin type with external white noise. The 
large equal-r equilibrium limit of the covariance function of the solution to this 
Langevin equation gives the one-dimensional two-point Schwinger function of the 
Euclidean field (j)a^j{t). 

The nonlocal Euclidean action of the massive scalar field which satisfies the 
(one-dimensional) fractional Klein-Gordon equation 

[(-A)" + A2]^0„^^(i)=O 

is given by 



(5.1) 



S[^an] = \ I 'f'cnit) [{-^r + X^V<l>c.nit)dt. 



Parisi-Wu quantization procedure requires (/)Q_^(t;r) satisfying the following sto- 
chastic differential equation 

d(j)a,^{t;T) _ 5S[(l)a.'f] 



(5.2) 



dr 



CKiT 



+ vit;r), 



2,.-,=<Pa,~,{t\T) 



where •q{t, r) is the external white noise defined by 

(r;(i; r)) = 0, (^^(t; r),7(t'; r')) = 25{t - t')5{T - r'). 

Eqs. ((5T|) and ((5?2|) give 



(5.3) 



^^^g^ = - [(-A)" + A^] -^ 0„,,(i) + n{t, r). 



The solution of (|5.3|) subjected to the initial condition (j)a^^{t; 0) = is 



(5.4) 



(t>a,jit;T) 



G{t-t';T^T')ri{t';T')dt'dT', 



where G(i; r) is the retarded Green's function given by 



G{t;T) 



(r) 



2tt 



exp 



(l-P" + A^) 



*duj. 



The large equal-r limit (t = t' —^ oo) of the covariance function gives 
(5.5) lim ((?!'a,7(ii,ri)(/)Q^(t2,T2)) 

Ti=r2^oo 

=2 lim // G{ti-t';T-T')G{t2-t';T-T')dT'dt' 



exp 



2a , \2\7^/ 



-2(|c.r + A2) 



e"^^*'-''^dT'duj 



^iu){ti—t2) 



2n J (|u;|2" + A2)' 



rduj, 



which is just the Euclidean two-point function of the fractional field (pa^^it). 

We next consider the Parisi-Wu quantization method at positive temperature 
T — 13^^ ■ We follow the Matsubara imaginary time formalism of finite temperature 
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field theory by requiring (/)^ ,^(i; r) to be periodic in the EucUdean time t with period 
p. That is, 



(5.6) 



€At + P;T) = 4,l(t;T). 



In addition, the white noise 77"^ (i; r) is assumed to satisfy the periodic condition in 
t, i.e. 77-^ (i + (3]t) — r]'^{t; r), such that 

(5.7) 

(77^(i;r))=0, (7?^(i;T)7y^(t';r')) = ^ E exp [zc^^t - t')] <5(r - r'), 

n— — 00 

with ujn = 2nn/p. The fractional operator [(— A)^" + A^] acting on (f>^ (t^r) is 
defined via Fourier series expansion (with respect to t) of cf)^ {t; r), i.e.. 



n=co r „^ 
[(-A)^" + AY0l.(^;r)^i E / <^..,(i' 



T)e-"^"''dt' 



n— — 00 



|2a I \ 2\ '^ iuJnt 



(KP- + A2) 



The retarded Green's function for the Langevin equation (|5.2|) satisfying the peri- 
odic conditions is given by 



(5.8) 



G^{t;r) 



ir) 



J2 ^^p 



IcoJ'- + X'yr 



Ji^rzt 



n— — oo 



The solution with initial condition (j>^ {t; 0) = is 



(5.9) 



€At;r) 



p r 



G'{t-t'-T-T')'n'{t'-T')dT'dt\ 



with covariance given by 
(5.10) 



^l.,^{tl.Ti)€,^{t2,T2)) ^ - Y. e 



iuJn[t\—t2) ^ 



n^ — oo 



L 



TxAT2 



exp 






exp 



(Kp" + A2)^(r2-r') 



dr' 



At the large equal-r limit, the covariance becomes 



(5.11) 



1 



Jirn^ (0^,^(il,Ti)0j_^(t2,T-2)) = - E 



3i(J„(tl-t2) 



ri— T2 — >oo 



/^„e^(KP"+^')" 



n— — 00 



which is the thermal two-point function for the Euclidean fractional Klein-Gordon 
field. When a = 7 = 1, this reduces to the ordinary two-point function of one- 
dimensional Euclidean scalar field at finite temperature [37] . We would also like to 
mention that when /? — > 00, the limit of the thermal two-point function (|5.1H) is 
the two point function (|5.5p . 

We remark that we can also consider the solution 0^^(i;T) to (|5.2p satisfying 
the initial condition <p'^At] —00) = 0, i.e., instead of having the field evolves from 
r = 0, we require it to evolve from r = —00. The solution 4''^^~f{t]T) is then given 
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by 

(5.12) (j)^.^(t;T)= I I G{t-t';T-T')rj{t';T')dt'dT'. 



It can be shown that (j)"^ {t;T) is a stationary field and in the large r limit, the 
field (j)a,j{t;T) approaches the field 0^^(t;T), i.e., 

lim 0a,^(i;r) ^cj)^ (t^r). 

The equal-r variance of 4>a^.y{t; r) is independent of r and is precisely the propagator 
(j5.5p . The same statement applies to (p^'^it; r) which is the solution of the periodic 
version of (|5.2p with boundary condition (jy^'^it; — oo) = 0. 

We would also like to remark that just as the stochastic process Xa^i3{t) is related 
to the field (jja.pit) in the sense that the covariance function of Xa^pit) coincides 
with the propagator of <j)a,p{t), we can define a periodic stochastic process X'^ ^^{t) 
whose covariance function is the propagator of the periodic field (/>^ (i) (|5.1ip . In 
fact, consider the solution of the fractional stochastic differential equation 

[(„Ar + A2]*Xj^,(t)= 77^(0, 
where rj^ {t) is the periodic white noise with period (3 and 

. oo 



n— — oo 



(r7^(t))=0, {^T{t)^T{t')) = - > e 



Using Fourier series, it is easy to check that the solution is given by 
and the covariance function is 



{xi^ml.it')) = 



[3 ^ (|c^„|2"+A2)^' 



n— — OO 



which coincides with the thermal two point function ()5.11|) . 

Now we proceed to calculate the partition function for the fractional oscillator 
process or fractional Euclidean field in one-dimension at finite temperature, hence 
its associated Casimir free energy. For this purpose, we employ the technique of 
zeta function regularization [38l [39l [40l [41]. Due to the fractional character of 
the scalar field under consideration, the derivation of Casimir free energy is more 
complicated as compared with ordinary scalar field. 

By definition, the Casimir free energy T of the fractional Klein-Gordon field 
(/)^ (i) which is kept at thermal equilibrium with temperature T — /3~^ is given by 

P 
where Z is the partition function defined by 

Z = JP^^ exp ("-i ^%^(i) [(-A)" + m^]^ ^^WdiV 
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Using zeta regularization techniques, we find that 

(5.13) ^=-^(C'(0)-C(0)logM'), 
where /^ is a normalization constant and C(s) is the zeta function 

oo 

(5.14) C(s) = m-27- + 2"^ {[an]^'' + m^y" , 

n=l 

with a = ^. The series in (|5.14p is divergent when s < l/{2aj). Therefore we 
need to find an analytic continuation of ^(s) to a neighbourhood of s = 0. For this 
purpose, we use standard techniques and write 

(5.15) 2y |[an]2"+m2r^^=-i- / t^^"! V e-*(['^"l"+™')dt. 

Since the obstacle for this integral to define an analytic function in s comes from 
the singularity at i = of the integrand, the asymptotic behavior of 



2E' 



-f:[an]" 



as ^ ^ becomes crucial here. Using the representation 

pc-\-ioo 

(5.16) 



1 pc-\-ioo 

mi J c—ioQ 



c e 



we have 



CO or 



c-\-ioo 



dwT{w)t-'"a-^°''"C,R{2aw), 



when c > \/{2q). Here C-r(s) is the Riemann zeta function. This gives the asymp- 
totic behavior 



(5.17) 5^, 



^r(^)™-.2E 



k\ 



t a CK(-2afc) 



n=l ^ ' k=0 

as t ^ (also as a ^ 0). Using the fact that Ci?(0) = —1/2, then with 



1/1 

if(i) = 2^e-*[""l'°--r' 

n=l 



a \2a 



i-^TQ-i + l 



eq. (|5.17p implies that K(t) ~ 0{t) as t — > 0. Now, we can continue the evaluation 
of the integral in (|5.15p : 



rj /'OC CG 

r(7.) Jo h 



r(7s; 

1 
^r(7s) Jo 
_ir(^)r(7s-5tj)^_, _ 



^7^ 



-r( — u-^fl-^-i 

a \2a 



r(7s) 



e-*™rfi+_l_/ t'''-'K{t)e''"''dt 



r(7s) 



a "m 2T"*+i _ 771-2-^" 



1 



r(7s) Jo 



t^''^K{t)e'''^ dt. 
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Since K{t) ~ 0{t) as t ^ 0, the second integral in the last line of the above equation 
defines an analytic function for s > — 1/7. Combining with the first term in (|5.14p . 
we find that an analytic continuation of ^(s) to Re s > — 1/7 is given by 



r(2b)r(7.-o^ 



2"^a-im-2^^+^ 



(5.18) as) ^, , 

To evaluate C,{0) and C'(0)j we observe that 



f<"'^K{t)e- 



dt. 



r(7s) Jo 

;hat 



is analytic for s > — 1/7. Therefore the only possible contribution to C(0) comes 
from the first term in (|5.18p when l/{2a) G N. Denote by A the set 

1 



A = <! — : ?i e N 
2u 



and let 



Then we find that 



i^a.A 



1, if a e A 



0, if a 1^ A 
C(0) = 2wc..A(-l)^a"^m = 



and 



C'(0) =2oj^^A{-l)^a~^miy U (^ + A~ ^(1) " log?™' 



27r 



7(l-u;a,A)^^pa-im- +7 / t ^K{t) 



2a 



dt. 



Here 4'{z) = ^'{z)/T{z) is the logarithmic derivative of the gamma function. Sub- 
stituting the above into (|5.13p gives the Casimir free energy 

^ if . ,N 1 ^ 1 



\og^?~-1[i:[ — + l] -i>{l)-\o; 



(3 



+ 7(l-w„,A)^^mi -7 / t-'^Kit) 



dt 



To study the asymptotic behavior of J-', we first use (|5.17p to obtain an asymptotic 
behavior of ^(s) when a —f 0: 



^(^) ^ r(2^)r(7^-^) 



aT^-ys) 



■a-'m-^-'' 



-U2^ 



(_i)fer(7. + fc) _2^ 



s-2k2ak 



fe=l 



kl 



r(7s) 



a'"''CRi-2ak). 



From this we can find the asymptotic behavior of C(0) and C'(0), which, when 
substituted into (|5.13p gives 



(-1)5^ 1 

(5.19) :F ^ui^^A^^ — TO" 



2tt 



log^t -^ 1 ^ I ^ + 1 



V-'(l) — logTTl" 



+ 7(1 — UJn, a) :^ 

'^ ■ ^2sin^ 

2a 



fc=l 



(-1)* 



\2afc 



m-^'^(27r)""'^Cfl(-2afc)/3 



-2Qfe-l 
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when /? ^ oo. In the special case a = 1, since ({—2k) = for all fc G N, (|5.19|) 
gives us 

(5.20) T^^. 

In fact, by applying the Jacobi inversion formula 

oo _. y — oo 



a V ^ 

n— 1 n— — oo 



we find that when a = 1, 



ar(7s) Jo „f:'^ 



n— — cxD 



aT{^s) aT{^s) ^-^ \amJ '^^^' 

^ ' ' ^ ' ' n— 1 



Here Kv{z) is the modified Bessel function of second kind. Together with Kx/2{z) 
\jTT/{2z)e^^ , one obtains 

(5.21) 



I3^n 2 /3 °^ '13 



n=l 



2 sinh 

2 



It is easy to see that the leading term when /3 ^ 1 agrees with (15.201) and the 
remainder terms decay exponentially. 

From (|5.19p . we note that at low temperature T ^ 1 (/? S> 1), the leading order 
term of the free energy J^ is of order T". When a ^ A, i.e. a is not the reciprocal 
of an even number, then the leading order is 

(5.22) ^^_2!^ + o(Ti+2"), 

2q 

whose sign depends on a. There is a dependence of J-^ on the normalization constant 
fj, when a G A. We shall renormalize the free energy J^ to get rid of this dependence 
later. 

To study the high temperature behavior of J-, we can use the expansion 

(5.23) exp(-tm2) = ^ i-^ 

1=0 

and find that 
(5.24) 

n=l ^ ^'^' 1=0 -^0 n=l 

2 ^{~iy 



I'm"' 



J2 ^m2'a-2a(^s+0r(^, + 0Cfl(2a(7s + 0)- 



r(7s)f- /! 



i=0 



The I = term in (|5.24p will contribute —1 to ((0), which is canceled by the 
contribution from the term m^'^''^ (|5.14p . Moreover, since C,r{s) is meromorphic on 
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with a simple pole at s = 1 with 
1 



CRis) 



s- 1 



^(l) + O is - I) ., 



as s ^ 1, 



we see that if a = l/(2j) for some j G N, then there is another nonzero contribution 
to C(0) arising from the I = j term in (|5.24p . Therefore, 



C(0) = 2w„,A(-l)^a~^m° 



and 



C'(0)=-27logm-2a7log— + 27 V LJ±m'^i a'^"^ Cni'^al) 

a ^-^ I 



2cc;„,A(-l)^a-^m^7(a|loga" + 2^(1)1 - V- ( ^ ) +^(1) 



This gives us 

(5.25) 



^ =^< 7logm2 + 2a7log/3 - 27 ^ 



^""(ll f"'^"" 



+ ^a,A(-l)='°-"^ 



log /i + ^7 log 



2V(1) 



^ 



2a 



V(l) 



In particular, in the high temperature limit T S> 1 (/3 ^ 1) 



(5.26) 



jT a7Tlogr + ^ logTO^ + 0(1), 



the Casimir energy is negative and the leading term — a7nogT depends linearly 
on a and 7. When a = 1, using 



2(20! 



32;, 



where -62; is the Bernoulli numbers defined by 



1 



1 1 



B 



1x2 



^ (27l) 



2n 2n~l. 
X 1 



and 



log 



l-e"^ 



2: Jo I e" 

we find that (|5.25p gives 



1 u 



(i?i 



X 1 



E 



D ^2n 

-D2n a; 



1 f °° 

■^ "25 1 "^^^S™' + 27 log/3 + 7^ 



=^<J21og[/3m]+21og 



1-e 



-,3m 



/3m 



2 2 ^^ (2n)! n ' 

n— 1 ^ ^ 

2 sinh 



(21)1 I 



f3m} = -log 



agreeing with (|5.2ip . 

One notices that when a e A, the free energy depends on the normalization 
constant /i. In order to remove this dependence, we need to renormalize the free 
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Figure 7: The renormalized free energy ^rcn/n^° as a function of I3m^'°' when a 
0.1,0.3,0.5,0.7,0.9 and7 = 1 



Contour plot of ttie function F /m as a function of p m and a 




Q r^ iatoffl^^Q itviaT^.tnJaiag;, rrvaT^lnA lSaK ,rrYtn^ 
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pm"" 



Figure 8: Contour plot of the renormalized free energy Trcn/m" as a function of p-m}'"' 
and a. Here 7=1. 



energy by adding a counterterm Tc to the free energy so that the renormahzed free 
energy J^ron is 

A reasonable way to determine the counterterm Tc is to require that in the limit 
/3 — !■ cxD and to -^- 0, T^cw -^ 0. Eq. (j5.19p gives us immediately 



Tc = — cJq,a 



27r 



-m° (log/i +7log?7i ) 



Note that adding the counterterm Tc to the free energy is equivalent to setting 
\x = m~^ when a e A. Therefore, from (|5.19p and (|5.25p . we obtain immediately 
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Figure 9: The renormalized free energy Trcn/rna as a function of (3m^'°' and a. Here 
7 = 1. 
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Figure 10: Left: This figure shows the ratio of the renormahzed free energy J-icn to 
— a7nogT when /3 = 1/T is small. Here 7=1 and m = 1. Right: This fig- 
ure shows the small T (large l3) behavior of the renormalized free energy .7-"rcn. Here 
Kc,-y = 2sin^/(7m^/°) and 7 = m = 1. 



that in the low temperature T = 1//3 ^ 1 limit, 
(5.27) 



(-1)^ L f , f I 

■?^rcn '^ - 7^Q,A 7; m" hA 7^ + 1 



^(i; 



2sm75— 

2a 






^-!^m-2'=(27r)2"'=C^(-2Qfc)/3-2"'=-i; 
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whereas in the high temperature T = 1//3 ^ 1 hmit, 
(5.28) 



•^rcn =^ S 7 log m^ + 2a7 log /3 - 27 ^ ■ 



(-1)' 2.^/3 



; V 27r 



2al 



^ Cfl(2a0 



+ 7'^a,A(-l)^° -TO° 



TT 



27r \ ^ . , , . \ f I 



«|log(^ +2^(l)l-V.(^-)+V.(l) 



Notice that (|5.27[) impUes that if a e A, then when T ^ 0, the leading order term 

of J^rcn is 



(5.29) .F,en ~ -r-^^rn^ U ( ^ + 1 ) - V^l) 1 + 0{T'"+'), 

whose sign depends on a. 

From ()5.22|) and ()5.29|) . we see that in the low temperature limit, the leading 
order term of the renormalized Casimir free energy still depends linearly in 7, but 
its dependence on a is highly nontrivial. 

Finally, we observe from (|5.27p and (|5.28|) that the combination T^cn/m^ de- 
pends on (3 and m in the combination I3m~ . The graphs and contour plots of 
J^ren/'rn^ as a function of a and /3ma are shown in Figures 7, 8 and 9. The high 
temperature and low temperature behaviors of J^^cn are shown in Figure 10. 

6. Concluding Remarks 

We have introduced a new Gaussian process called the fractional oscillator pro- 
cess with two indices, which is obtained as the solution to a stochastic differential 
equation with two fractional orders. Some basic properties of this process can be 
obtained based on the asymptotic properties of its covariance despite its complex 
nature. The possible generalizations of fluctuation-dissipation theorem and Ein- 
stein relation to this fractional process are considered. The main advantage of the 
oscillator process parametrized by two indices over the process with single index is 
that the former has its fractal (or local self-similarity) and the short range depen- 
dence separately characterized by the two indices, while the later has both these 
properties determined by a single index. Such a process may provide a more flexible 
model for applications in phenomena with short memory. 

By analogy regarding the fractional oscillator process as the Euclidean fractional 
scalar Klein-Gordon field in one dimension, we carry out the stochastic quantization 
of such a field with a nonlocal action. The Casimir energy associated with the 
fractional Klein-Gordon field at finite positive temperature was calculated using 
thermal zeta function regularization technique. The expression for the free energy 
has a rather complicated form. We thus consider the low and high temperature 
limits for the free energy. Graphical representations of these asymptotic behaviors 
of the free energy are given. 

Extension of our results can be carried out to give the n-dimensional Euclidean 
fractional Klein-Gordon field. However, the derivation of the asymptotic properties 
for the covariance function and the Casimir energy will be more complicated. The 
sign dependency of the free energy on a has important physical implication when 
the fractional quantum field under consideration is confined between parallal plates 
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or cavities as the sign of Casimir energy will determine whether the associated 
Casimir force is attractive or repulsive. 

Another direction of generalization is that of fractional oscillator process of vari- 
able order, with a and 7 being extended to time-dependent a{t) and j{t). We 
expect the variable short range dependence property remains valid, and the result 
for fractal dimension holds only locally. However, the Casimir energy calculation 
will require new mathematical techniques and approximations. Such a generaliza- 
tion may find applications for complex systems where the physical phenomena can 
have variable short memory and the fractal dimension varies with time or position. 
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Appendix A. Mathematical details 

1. We show that if the position process Ya^.y{t) is related to the velocity process 
Xa,-yit) by 



Xa,^{t) 



''''' - dt ■' 



then the long time behavior of the covariance function Ca.-f{T) of Xa.-y{t) does not 
show up in the leading term in the long time asymptotic expression of the variance 
[i^a, 7(^)1 )• Its effect only appears in the second leading term. In fact, from (14. 7p . 
we have 



{[yo.,^{tf) ='2 f (t- T)C^^^{T)dT = 2 



Ca,'y{T)dT 



t-2 TCa,-y{T)dT. 







Eq. (|3.18p gives us the long-time behavior of Ca.^ir) as Ca^-yir) ~ Ait~'^°'^^ + 
0(r^'*"~^) for some constant Ai, which implies that the integrals /q Ca,^{T)dT is 
convergent; and 

(A.1) 

pt POO poo 

/ C„,^(T)dT = / C^.^iT)dT - / Ca,,^iT)dT ^ A^ + ^3*"^" + 0(^-4") 

Jo Jo Jt 

as i — !■ 00. On the other hand, the integral L TCa^^{T)dT is convergent if and only 
if a > 1/2. In this case 

(A.2) TCc,,y{T)dT= TCa,j{T)dT- rCa,7(r)dr~ A4 + 0(r^"+^) 

Jo Jo ' Jt ' 

as t — + 00. If a < 1/2, then 

(A.3) ( TC{T)dT - Agt^-^a _^ Q(^max{0,l-4a}) 

Jo 
as t ^ 00. In the borderline case a — 1/2, 

(A.4) / TC(T)dT- A log t -I- 0(1). 

^0 
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Therefore, as i ^- c«, if a > 1/2, then 

If a = 1/2, then 

([Y^,y{t)f) ^ [2A2]t - [2Ae] logt + 0(1). 
Finally, if a < 1/2, then 

'[Y^^^{t)f) ^ [2A2]t + [2A3 - 2A]tl-2a + O(imax{0,l-4a})^ 



These show that the leading term of ( [Vq .^(t)] \ is of order t, independent of a; 
and the second leading term is of order imax{o,i-2a} logi, which depends on a. 



2. We show that if the position process Ya,-y{t) is related to the velocity process 
Xa,-yit) by 



Xa,-f{t) — O^t ^Q,7(0: 



2<^<2' 



and QD}~^Yo,,^{t) 
(A.5) 



t=o 



for j = 1 if X < 1 and j = 1 , 2 if x > 1 , then 



[r„^^(t)]n = 2i?A-2"7 



tix- 



{2x - \)V[x? 



+ O /'imax{0,2x-2,2x-2a-l} j^g A 



From (|4.1ip . we have 



[^".7W]')=w^/ / (t-Si)X-l(i-S2)>^-l(X„,^(si)X„^^(s2))dsldS2 



nt pt 

Using some calculus, this gives 



^ / /(t-Si)X-l(i- 52)^-1^,^(151 -S2|)rfsldS2. 

X) Jo Jo 



(t-SiY ^{t~S2)^ '^Ca,^{s2~ Si)dsids2 



nx)' Jo Jo 

nt nS 



nxY Jo Jo 



{t + T- s)^-\t - sY-^Ca,,^{T)dTds 



nx)'- 



■ / / {t + T-sY-\t-sY-^ds 



Ca,-y{T)dT 



,x-i 



[u — t)^ ^du 



Ca,-y(j)dT. 



nxY Jo 

The case x = 1 ti^s been considered above. Now we consider the cases x G (1/2, 1) 
and X G (1, 3/2) separately. If x ^ (1/2, 1), using 



(u - t)x-i = u>^-i - (ux-i - (li - t)>^-i) 



we have 

M-^" (w — t)-^~ du : 



2X- 



- (t^X-l _ ^2x-l) _ /■ yX-1 (j,X-l _ (y _ ^)X-1) rf^_ 
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Therefore, if x e (1/2,1), 



CJ.2X-1 /■* 1 /■* 



liX-i fuX-i 



Ca.-y{T)dT 



On the other hand, if x G (1, 3/2), using 

{u - t)^-^ = u'^-^ - (x - l)u^-^T - [uX-1 _ (;^ - l)uX-V - (li - t)'^-i] 
we have 



i^~\u-T)^''^du 



Therefore, if x e (1,3/2) 



t2x- 



2x-l 2 



-t^^ 



-2 _ 3-2x ^2y-i 

2(2x-l) 



^X-l [uX-l _ (;^ - l)wX-V - (u - t)^-^] du. 



2 / +2X-1 rt f2x~2 ft 

3-2x /■* 2^-1 



2(2x-l)yo 



'-^^"^C„^^(T)dT 



I / U^-l [u'^-l - (X - 1)UX"V - (U - t)^-^] du\ Ca,-yiT)dT 



k"'r 



The large-i behaviors of L Ca,-f{T)dT and L TCa,-f{T)dT have been studied and 
given by jXTj) . ([X2)) . (fO]) and (jXI)) . For the term /J r2x-iCa,^(T)dT, since 
Ca.7(T)(T) ~ Air-2"-i + 0(t-4"-1) as t -> oo, we find that if 2x - 1 < 2a, 



T2^"^Ca,7(T)dT = / T^^-^Cc,,^{T)dT~ T'^^-^Cc.,^{T)d 



However, if 2x — 1 > 2q;, 



=Bx + 0(^"^"+^^"^). 






and if 2x — 1 = 2a, then 



/ T'^^-^C^^^{T)dT^B^\0gt + 0{l). 

Jo 



Now for the term 



,X-1 /'i/X-l 



(uX-i-(u-T)>^-i)du 



CQ,7(T)dr, 
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By making a change of variable u '—^ ut, we have 



,x-i (o,x-i 



{u — t)-^ ) du 



Ca.-f{T)dT 



,2x-l 



,x-i A,x-i 



{u-l)^-^)du 



Ca.f{T)dT 



x-1 /^x-i 



(w-l)X-l)/ T^X-^Ca,^^{T)dTdu. 

Jo 



Using the fact that Ca,-y{T) > 0, we find that this term is bounded above by 

CO pt 

yX-1 (yX-1 _ (y _ i)x-i) du / r^x-i^^ ,^(^)(^)^^_ 
1 -/o 

Notice that x ^ (V^, 1) imphes that the first integral is convergent. Similarly, we 
have for xe (1,3/2), 

ft r rt > 

uX-1 ^yX-l _ (^ _ l)u'^-V - (m - r)'^-l] du \ Ca.,-y{T)dT 



< 



,x-i r„x-i 



[m^~^ - (X - 1)"''"^ - (" - 1)''"^] du / T2^-lCa,^(T)d- 



'1 JO 

and X € (1,3/2) guarantees the convergence of the first integral. Gathering the 
results, we find that when t — > cx). 



(^[Y^^^{t)f) ^ At^^-' +0 {mi,x{t^^-^''-\t^^-\t^^-Hogt,t°,t°\ogt}) 



where 



A 



(2x - i)r(x)2 Jo 



Ca,'y{T)dT 



2SA-2"T 

(2x-i)r(x)2' 



and (JA.SP follows. 
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